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Abstract 

Quantum extensions of the Gallavotti-Cohen fluctuation theorem (FT) for the entropy production have been 
discussed by several authors. There is a practical gap between microscopic forms of FT and mesoscopic (i.e. not 
purely Hamiltonian) forms for open systems. In a microscopic setup, it is easy to state and to prove FT. In a 
mesoscopic setup, it is difficult to identify fluctuations of the entropy production. (This difficulty is absent in the 
classical case.) We discuss a particular mesoscopic model: a Lindblad master equation, in which we state FT and, 
more importantly, connect it rigorously with the underlying microscopic FT. We also remark that FT is satisfied 
by the Lesovik-Levitov formula for statistics of charge transport. 



1. Introduction 

A goal of nonequilibrium statistical mechanics is to make good, yet simple, models of complicated phenomena 
out of equilibrium and to analyze them. Such models should answer correctly questions like: "What is the structure 
of an object's state that is kept between two heat reservoirs at different temperatures? How does it look like locally 
and how does it fluctuate?" or "When and how does a liquid flow become turbulent?" . In equilibrium, we have a 
good and simple guiding rule for making dynamical models: detailed balance. Surely, any nonequilibrium dynamics 
should break detailed balance, but can one formulate a general rule as to how to fix this detailed balance breaking? 

It is here that the celebrated Gallavotti-Cohen fluctuation theorem (FT) (see [18,19,20]) for the entropy pro- 
duction comes into the picture. It turns out that there is such a general rule and it reads (for every path a of the 
dynamics, with corresponding probability Prob[<r]) 

Problcr] 

log : — - — Entropy production along a + o(t) (1) 

Prob [6* <r J 

where 9a is the time reversed path a. The FT itself 

Prob [Entropy production up to time t equals st] , 

log — — ■ ; ; -r; = St + 0(tj (2) 

Prob [Entropy production up to time t equals — st \ 

is an immediate consequence of the above rule (1). The LHS are expected to be linear in time t, the o(t) term 
vanishes after dividing by t and taking t f oo. See [29] for a discussion focusing on the universality of (1). 
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For our purposes, it is good to emphasize that there are two reasons to appreciate (1,2), even if those reasons 
are obviously intimately related: 

(I) Formulas (1) and (2) can be derived on a microscopic level, starting with Hamiltonian dynamics, using its 

time-reversibility and identifying entropy production as the log of phase space volume. This was outlined 

in [30]. 

(II) The rules (1) and (2) do apply to a lot of effective (not purely Hamiltonian) or mesoscopic models. It is 
actually in effective models that the FT was observed and that it is studied most often. In [18,19,20], one 
treats thermostatted systems, while in e.g. [25,26,11], as well as in most later publications, one deals with 
stochastic models. 

Note that (II) does morally, but not exactly follow from (I) since mesoscopic or effective models are merely 
good approximations of the Hamiltonian microscopic dynamics. 

This article is concerned with the question how rules (1) and (2) extend to quantum mechanics, On the mi- 
croscopic level (cfr. (I) above), there is an obvious problem with deriving (1); quantum mechanics does not have 
the concept of paths with associated pathwise entropy production. However, (2) still makes sense and in fact it 
remains true, as was pointed out in [22,24,38,35] and as we review in a streamlined version in Section 2. 

As to the status of quantum FT in effective or mesoscopic models (cfr. (II)), the majority of research has 
focused on master equations for small systems (see recently [28]). A very explicit discussion can be found in [36], 
a simplified version of which is presented in Section 3.3. The ideas developed there are also implicit in [35] and 
[33]. 

The fact that not more effective or mesoscopic models have been investigated is perhaps not too surprising. With 
the possible exception of the master equations discussed in Section 3, there are not much Markovian models of 
open quantum systems which are healthy in all respects. One frequent pathology is the appearance of nonpositive 
density matrices. A review of this problem and a list of possible approaches is given in [12]. We should mention 
that these problems are not inherent to nonequilibrium models. Already the notion of detailed balance is somehow 
problematic for Markovian quantum models since the existing definitions [2,3] are so restrictive that the master 
equation in Section 3 is essentially the sole example. 

Moreover, even if one has a trustworthy model, there still remains to identify the fluctuating entropy production. 
It should be clear from Section 3 and Remark 1 in Section 3.3 that this is not trivial. 

We remark hence that in quantum mechanics, there is a practical gap between (I) and (II). This might be 
disconcerting, but possibly it also offers opportunities: perhaps one can improve upon some models of open 
quantum systems by imposing FT. 

In the case of master equations which we treat in this article, not only (2) holds but there is also a convenient 
framework to state (1). This framework is the formalism of unraveling of master equations. Our results imply that 
these unravelings are rigorously linked to fluctuations of currents. This is discussed extensively in [16]. 

1.1. Brief summary 

In Section 2, we review and sharpen the Hamiltonian approach to the FT. In Section 2.3, we give an example of 
a known counting formula in mesoscopic physics which satisfies FT. We introduce the master equation describing 
a small system in contact with different heat baths in Section 3 and we develop the necessary tools (unravelings) 
to identify the entropy production. Finally, in Section 4, we connect the master equation with the Hamiltonian ap- 
proach and we show that fluctuations of the entropy production in the Hamiltonian model converge to fluctuations 
calculated by means of unravelings. The importance of this result is: 

1) It shows that the FT is valid in the thermodynamic limit (at least in the weak coupling limit). 

2) It rigorously justifies the technique of unravelings, which is quite popular in quantum optics. 



2. The Quantum fluctuation Theorem 

To fix thoughts, we introduce a model and corresponding notation. Imagine several heat reservoirs indexed by 
k 6 K and modeled by Hilbert spaces Hk and self-adjoint Hamiltonians Hk on Hk- Each reservoir is in thermal 
equilibrium at inverse temperature /3k- One connects all heat reservoirs with a small system S, modeled by a 
finite-dimensional Hilbert space Hs and a self-adjoint Hamiltonian Hs- The connection between S and reservoir k 
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is through a self-adjoint coupling term Hg_ k on Hs ® fik- The coupling is controlled by a factor A and the total 
Hamiltonian is 

fc fc 

The dynamics of the full system is given by the unitary group Ut = e~' tHx on the Hilbert space Tt := 
T~ts ® [®kHk\- The initial state po is a product state on TL 



pa = po 1 



®Pk,p k 

k 



(4) 



where the states pk,p k are equilibrium states on TCk with respect to the Hamiltonians Hk at inverse temperatures 
j3k and pio is an arbitrary density matrix on Tis- 

In the thermodynamical limit, the states pk,p k are no longer density matrices and one needs the machinery of 
operator algebras, see e.g. [7,14], to define these states properly. We will completely ignore this problem (since its 
solution belongs to standard knowledge in mathematical physics), and we avoid writing ill-defined expressions by 
using the notation 

Pk[A] := Tr[pfcA] with A an operator and pk a density matrix on Hk (5) 

The LHS can still be used when pk is no longer a density matrix, but instead a state in the sense of operator 
algebraj^l Note finally that models as (3) have received a lot of attention lately in mathematical physics. E.g. in 
[21] one proves the Green-Kubo relations and Onsager reciprocity in the spin-fermion model at finite A > 0. 

2.1. What do we mean by quantum fluctuations? 



It is not a priori clear what one means by fluctuations of heat currents, or fluctuations of the entropy production. 
(See [5] for an elaboration on that question.) 

One encounters at least three approaches (see however [5] for a different view), 

1) Measure the energies Hk in the beginning and at the end of the experiments and make statistics of the 
difference of both measurements. This approach yields the FT [35,22,24,37]. 

2) Calculate the fluctuations of the operator UtHkU-t ~ Hk- In [32], a deviation from the fluctuation theorem 
was established for the fluctuations of a related operator. 

3) In [31], one studies the fluctuations of the relative modular operator - a C*-algebraic concept. These fluctu- 
ations satisfy the FT, but the meaning of this relative modular operator is a priori unclear to us. 

The advantage of the first approach is that it is more relevant in an experimental setup. A possible argu- 
ment for the second would be that it seems better suited for taking the thermodynamic limit, since it avoids 
the measurements. However, by slightly reformulating the first approach, one can also meaningfully study the 
thermodynamical limit (see also [23]), as we point out now. 

Assume for simplicity that the Hk have discrete spectrum, indicating that we have not taken the thermodynamic 
limit and let x 6 X label a complete set of eigenvectors \x > of (Hk)keK with eigenvalues (Hk)keK(x). Let the 
projectors P x stand for 1 <g) |a;)(a;| where 1 stands for identity on Tis- The probability of measuring an entropy 
production r is 

P*,A(r) = ]T po [p x Ut t P y U t x P x ] (6) 

The idea behind the formula is clear: measure (thereby projecting on the eigenstates x), then switch on the time 
evolution Ui', finally measure again (projecting on the eigenstates y). For later convenience, instead of studying 
the probability distribution P tj ^(r), we focus on the Laplace transform. 

F(K,t,\, P0 ) := JdF tA (r)e- Kr (7) 



2 In fact, the situation is even more complicated. For fermions, one can indeed find a subalgebra of B(7ik) -The Weyl-C*-algebra- on 
which the state Pk,fi k can be appropriately defined. For bosons, one has to invoke W^-algebra's and the Araki-Woods representation 
to give rigorous meaning to the model. 
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which, in the case when W t ,\(r) is given by (6), reads 

F( K ,t,\,p ):= Yl Po[P*Ut t P v U t x P x ] e - K ^ MHM - H " l ' )) 



x, y ex 



We now use that the initial state po is diagonal in the basis \x > to rewrite (8) into 



F(n,t,\,p ) = po 



(8) 



(9) 



This last expression is perfectly suited for taking the thermodynamic limit, since it remains well-defined when 
the operators H k have continuous spectrum and the state po is not given by a density matrix. (Which doesnot 
mean that the expression (9) is necessarily finite. After all, the physics could be such that the Laplace transform 
is infinite for some n.) 



2.2. Derivation of the quantum fluctuation theorem 

If one defines the fluctuations through (9), FT is easily derived. Assume that the initial state po is given by 



d 



-0k H k 



Z k := Tre 



d := dimHs, 



(10) 



hence the initial state on the small system is the trace state, or infinite temperature state. Assume further that 
the dynamics is reversible, i.e. there is an anti-unitary time-reversal involution T which commutes with both the 
free and interacting Hamiltonians 



[T,H k 







[T, H x ] = TT = 1 TUtT = U_ 



Abbreviating W = ~}2 k f3 k H k , using the above relations and cyclicity of the trace, one calculates 
F(n,t,\,po)=Ti [e~ w TTU t x TTe~ KW TTU* t TTe KW j (df^) -1 



= Tr\e- vv U t A e 



te (i-K)wj (dJT^r 1 =F(l-M,A,Po) 



(11) 

(12) 
(13) 



which is equivalent to (2) without the o(£)-correction term. The above derivation relies on a particular initial state 
for the small system. It is known as a transient FT. Obviously, for large t, one expects the initial state of the 
small system to be irrelevant and (2) follows under some additional ergodicity assumptions. 



2.3. Example: counting statistics for charge transport 



As argued in Section 1, it is interesting to investigate to what extent the fluctuation symmetry is actually present 
in known expressions for current statistics. We briefly discuss the Lesovik-Levitov formula [27] for electronic 
transport between leads. We base our treatment on a neat version of the formula presented in [23] wherein the 
electrons do not interact and the formula (16) can be derived without any approximation. It is also assumed that 
all Hilbert spaces are finite-dimensional (see [6] for a discussion of the thermodynamical limit). 

Let f) be the Hilbert space of one electron. The presence of two leads is made explicit by splitting () = f)i © ()2 
where the spaces f)i, f}2 contain all the states in respectively lead 1, 2. The projectors on f)i, 1)2 are denoted 1^ , 1(, 2 
and the decoupled one-particle Hamiltonian is hi © fe. Let the inverse temperature be j3 and let the leads have 
respective chemical potentials pi, p2- Let n be the occupation number operator (dictated by Fermi-Dirac statistics) 

n = (V (hl - Mll "i> e P( h 2-v2U 2 ) +1^ 1 (14) 

The contact between the leads is modeled by a scattering process with unitary scattering matrix S. We assume, 
that the scattering process conserves total energy and that both scattering process and decoupled evolution are 
reversible 

[S,hi®h 2 ]=0 Thi, 2 T = hi,2 TST = S* (15) 
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for some anti-unitary operator T modeling time-reversal. 

At time t = and t = r, one measures the total charge in, say, lead 1 (For simplicity, we set the charge to 
unity, so that 'charge' just means 'number of particles'). Let F(k) be the laplace transform of the distribution of 
transported charge F(k) = J^J^-oo P(?) e ~ Kq where p(q) is the probability to find at time t — r an excess charge 
q, as compared to t — 0. 

The Lesovik-Levitov formula reads (see [23], where the formula is given for the Fourier transform, hence for 

k e iffi): 

F(k) = det (1 + n(S* e~ Kl »i Se Kl 'n -1)) (16) 

The entropy production in this scenario is just — /3(/m — ^2)9 and hence one expects that FT implies the following 
symmetry (cfr. (12-13)) 

- v*)) = F(-(l - - a* 2 )) (17) 

That relation is easily verified by starting from (16), using (15) and the fact that det(AB) — det(A)det(£>) for 
matrices A, B (Recall that our one-fermion space is finite-dimensional) . 



3. A good effective model: the master equation 

3.1. Construction of the master equation 

Having in mind the Hamiltonian model introduced in Section 2, we can construct a master equation that 
approximates the evolution of a general density matrix /10 on the small system Tis- This master equation describes 
essentialhPI two effects: 

1) Decoherence: non-diagonal elements in the eigenbasis of Hs vanish 

2) The fluctuations of the diagonal elements, which evolve independently of the off-diagonal elements. 

In view of this behavior, we choose to restrict attention to diagonal elements. This allows for a less technical 
presentation, as we can replace the Lindblad master equation by a classical Markov process. 

Let e,e',... be the eigenvalues of the system Hamiltonian Hs and let P e ,P e /,... be the projectors on the 
corresponding eigenspaces which we assume one-dimensional. 

Define transition probabilities for each seperate k £ K as 



Pk{e,e 



f At &V-'\l®p h ,p h ) [e iHht P e ,Hl_ h P e& - iH ^ P e H S - k P e '} (18) 



where 1 is the (unnormalized) density matrix 1 on Tis- The total transition probabilities are hence given as 

P(e,e') = ^p fe (e,e') (19) 
fe 

The information about the temperatures 0k is encoded into the RHS of (18) through the state pt, and in fact one 
can show 

Reservoir k is at thermal equilibrium at flu =>- pfc(e, e')/pk(e, e) = e' 3fc ' e_e ' (20) 

As said, the density matrix of the system is diagonal. We denote the time- dependent diagonal elements by p, £ (t). 
The (Pauli) master equation now reads 

= E WP( £ '' e ) " ^(iMe, e')) (21) 



That is, if wc assume Hs to be nondegencrate. If it is degenerate, the picture is slightly more complicated but the results still 
hold, sec [36] . 
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3.2. Derivation of the master equation 



The only information missing is in what sense this master equation (21) approximates the true dynamics 
generated by (3). Physically, one needs to perform the Born-Markov approximation and the rotating wave ap- 
proximation, on which we will not comment here (see instead [4,8]). Mathematically, one can make the convergence 
to the master equation precise in the so-called weak coupling limit: Let the initial state po = po ® [®kpk,p k ] be as 
in (4). Let p x (t) be the density matrix obtained by reducing the full state of the interacting system at time t, i.e. 

p x (t) = Tr R [[/ t A Mo ® [<8>kPk,p„] U- t ] (22) 

where TrR is a conditional expectation. (If p is a state on Hs <8> Wa, then TrR,p is a state on Hs such that for 
any operator A on Hs: p(A) = (Tr np)(A). If p is a density matrix, then TrR is just the partial trace over the 
reservoirs.) The statement, rigorously proven in [13] under mild assumptions (about which we will be more explicit 
in Section 4), is 

P e ^ x (\- 2 t)P e ^-> o p e (t) (23) 

As is clear from the scaling and the expression (18), the master equation emerges in second-order perturbation 
theory. Hence we require that TrR[i7<S_ fc ] = such that the first order term perturbation term vanishes (which 
can always be achieved by redefining Hs-) 



3.3. Guessing path-wise entropy production 

Given the very suggestive form of the transition probabilities (18), one is inclined to supplement the master 
equation in the following way: 

The paths a of the markov process defined by formula (21) are sequences of energy levels e, with transition 
times ti (jump from ej_i to e{) with Pg the path probabilities up to time t 

o = (eo; ti, ei; t 2 , e 2 ; . . . ; t„, e„) P* [a) = e^"* 1 p ( e0j ei ) e - A «i («2-«0 . . . p ( ej -i,e j ) e~ x °j (24) 

where j is the highest index such that tj < t and A e = J^ e , p(e, e') are the escape rates. The subscript eo in ¥ t eo 
indicates the the process has been started from eo. 

Clearly, to speak about currents and entropy production, we need more information; we need to know which 
bath k triggered the transition. Hence, we need a process with paths a and probabilities P* (<7); 

a = (e ; ti, fci, a; t 2 , k 2 , e 2 ; . . . ; t„, k„, e„) ¥ eo (d) = e" Ae o*i Pfci ( e0j ei ) e -^i('2-ti) ( ej _i, e^) e" Ae i ( *-*^ (25) 

We simply define a new process corresponding to the pathspace measure (25) . Remark that obviously, 

]TPl (a)=P*» (26) 

a — >a 

where the sum is over all a which reduce to a by omitting all ki . 

One can just as well calculate the matrix elements p e {t) from the path probabilites P as from P: 

where 

K(0) = E ( )^0 . P M(0) = E ^0 (0)Pe„ 

and with n(a),n(a) being the number of jumps in the paths a, a. Hence, we have extended our pathspace proba- 
bility measure P* into a pathspace probability measure P*. In the theory of open quantum systems, one calls such 
an extension an unraveling and the paths a are called quantum trajectories, see [9,8] for an overview and [10] for 
similar reasoning. 



(28) 
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Having the pathspace probability measure P* at our disposal, we can do all the manipulations which are familiar 
from classical markov processes. The function on paths a 

n 

w*(er) = -X)/3fci(ei-ei-i) ( 29 ) 

which we identify as the path dependent entropy production, is equal to the ratio of probabilities of forward and 
backward paths; 

lu'(cr) = log with 9a = (e„;t - t„, k n , e„_i; . . . ;t - t 2 , k 2 , e 2 ; t ~ ti, ki, e ) (30) 

dP w (6a) 

as follows from (20). From (30), we retrieve the FT (2), or its Laplace-transformed version 

log J dP^e-™'^ =log |dP^(5)e-( 1 -^ t (*)+o(t) (31) 

Obviously, one can push on and prove the Green-Kubo relations, Onsager reciprocity and strict positivity of the 
entropy production from (31). This is discussed at length in [36]. 

At this point, one should note that this unraveling and the construction of the pathspace measure P* was a 
product of our intuition. Hence, it not a priori clear whether the fluctuations of entropy production which we can 
calculate through formula (29) and the measure P coincide with the entropy fluctuations in the original model. 
This problem motivates the following section. 

Remark 1 However intuitive the above reasoning, one should keep in mind that we have used more information 
than is contained in the bare master eguation (21). In particular, we used (19). To put things even sharper: given 
a master eguation without any additional information, one could associate to it baths in different ways, thereby 
obtaining completely different expressions for the entropy production. 



4. Linking the microscopic fluctuations with the quantum trajectory fluctuations 

The idea of this section is to justify the guesses in Section 3.3. We do this by connecting fluctuations of the 
entropy production in the Hamiltonian model, as in expression (9) to the fluctuations calculated by the probability 
measure P* . This connection is made sharp in the regime in which the master equation is derived, i.e. in the weak 
coupling limit. To be as concrete as possible, we make explicit the model which was outlined in the beginning of 
Section 2. 

Let for each k, 7i k be a fermionicQ Fock space, i.e. 

Hk ■= C L 2 (R d ) © L 2 (R d ) ® a L 2 (R d ) © . . . (32) 

where the © a denotes the anti-symmetrized tensor product, which accounts for the fermion statistics. Let the 
Hamiltonian Hk describe free fermions with one-particle dispersion relations h k (q) = h k (\q\) where q € R d is the 
momentum, i.e. 

H k = J dqh k (\q\)aUq)a k (q) (33) 

with a k (q), a k (q) creators/annihilators of fermions in reservoir k with momentum q, that is 

{a k (q),al,(q')} = 5 k , k ,6(q - q) (34) 
Let now the interaction H\_ k be given as 

dqf k (q)D k <g> 4(g) + f k (q)D* k ® a k (q) (35) 



One could just as well take bosons here, but then formulas like (37) need more elaboration to get rigorous meaning, see also the 
footnote in the introduction to Section 2. 
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for some operators D k on Hs and coupling functions f k 6 L 2 (R d ). The thermal states pk,/3 k on Hk are defined as 
quasi-free states determined by their two-point correlation function (determined by the Fermi-Dirac distribution) 



Pk,p k [al{g')a,k{g)] = (g, (l + e Pk k ) g) 



(36) 



We recall the definition of F(k, t, A, po) in (9) as the Laplace transform of the distribution of entropy production; 



F(k, t, A, po) = Po 



: T, k /3kHk iX A e K £<= /3fc 



(37) 



where again, as in (4), po = po ® [®fcPfc,/3 fe ] with po an arbitrary density matrix on Hs- 
We can now state the main message of this paper. 



The generating function F(n,t, A, po) (37) /rom </ie microscopic Hamiltonian model con- 
verges in the weak coupling limit to the generating function defined by unravelings, 

F(n,t,IM>) = J dre- Kr fl [w* = r] 

R 

where Pj, is ifte measure constructed in Section 3.3. That is, under some conditions (see 
below) and for allt>0 and initial density matrix po on Hs, 

F(k, A~ 2 t, A, po) — > F(k, t, po) 



Rephrasing the fluctuation relation within the Hamiltonian model and taking t f co one gets : 

lim lim t~ x log F(n, \~ 2 t, A, po) = lim lim log F(l — n, \~ 2 t, A, po) (38) 

t/-|-oo A\0 f/ico A\(] 

which is a mathematical statement about large deviation generating functions. 

The limit A \ can be omitted altogether from the above statement since the FT holds without small-coupling 
approximation (see [34] for a proof of existence of the limit t f +oo in that case). 

These statements are essentially a rephrasing of technical results in [15], and they are more extensively com- 
mented upon in [16]. The results are proven by extending the convergence in the weak coupling limit to the full 
system (including reservoirs) . One proves that the full Hamiltonian evolution converges to a quantum stochastic 
evolution. Similar results (though in a weaker form, not permitting the conclusions of this paper) have been proven 
in [17] and [1]. 

To make the discussion complete, we give an example of a set of assumptions that is sufficient to prove the 
above statements. 

Define the set of Bohr frequencies T := {e — e' |e, e' € spffs} and assume that 

1) The dispersion functiosn hk : R + i— > R + are monontone increasing, and continuously differentiable with 
nonzero derivative in /i^" 1 (^ r ). The coupling functions fk are also continuous in ft^ 1 (J r ). Write hk for the 
self-adjoint operator associated to the dispersion function h k by (h k g)(q) = h k (\q\)g(q) for g G L 2 (R d ). 

2) For 7 = 0, k, the functions 

f p fc a, \a(e- [thk e -~< f3khk f)a*(f)] 
{PKP k [a(e- uhk e^ hk f)a(f)] 

are bounded and L 1 -integrable 



5. Conclusion 



We investigated the question of generalization of FT to quantum systems. We remarked in Section 2.2 that in 
a microscopic (=Hamiltonian) framework, one can unambiguously state FT. In practice, it is more common to 



8 



start from a mesoscopic (=effective, nonhamiltonian) description of an open quantum system. In Section 3.3, we 
discussed a very straightforward and successful approach for the dynamics given by a specific master equation; 
one can guess what are fluctuations of entropy production and one can state FT. Finally, in Section 4, we state 
the link between fluctuations of the entropy production, defined in the microscopic model in Section 2 and the 
fluctuations which are the result of the guess-work in Section 3, which hence turns out to be good guess-work. 
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